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Abstract. The millisecond proto-magnetar model for the central engine of long-duration gamma- 
ray bursts is briefly reviewed. Limitations and uncertainties in the model are highlighted. A short 
discussion of the maximum energy, maximum duration, radiative efficiency, jet formation mecha- 
nism, late-time energy injection, and (non-)association with supernovae of millisecond magnetar- 
powered GRBs is provided. 
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INTRODUCTION 

Models for the central engine of long-duration gamma ray bursts (GRBs) are highly 
constrained by the character of the prompt emission and the afterglow, and — at least in 
some cases (and perhaps nearly all yj]) — the fact of an associated supernova (SN). A 
successful model must at minimum satisfy several criteria. It must generate a collimated 
outflow with high Lorentz factor y TC ^ 10 2 — 10 3 and high kinetic luminosity {E ~ 
10 50 - 10 51 ergs s ') on a timescale of ~ 10 — 100 seconds. Additionally, inferences 
from X-ray observations by Swift indicate that there may be late-time energy injection 
by the central engine on timescales of hours to days Q2Q. 

There is little diversity among models for the central engine and essentially all can be 
simply classed as a rotating compact object that drives an asymmetric relativistic out- 
flow. The two leading models are the "collapsar" and the "millisecond proto-magnetar." 
The former, as described in ^ and developed in proposes that GRBs arise from 

the collapse of rapidly rotating Type-Ibc progenitors. A black hole forms with an ac- 
companying accretion disk and drives a collimated relativistic outflow along the axis of 
rotation via either neutrino heating or magnetic stresses. Ref [6] have recently pointed 
out that rapid mass loss during rotating collapse might prevent the formation of a central 
black hole, at least for some time. 

In contrast to the collapsar mechanism, the "millisecond proto-magnetar" model 
posits a newly formed rapidly rotating neutron star (spin period P ~ 1 ms) with sur- 
face magnetic field of magnetar strength (B ^ 10 15 G), cooling via neutrino radia- 
tion, and driving a neutrino-heated magneto-centrifugal wind 10, 0, @]. Millisecond 



proto-magnetars might be produced by rotating Type-II and Type-Ibc progenitors, the 
accretion-induced collapse of a white dwarf, the merger of two white dwarfs, and/or 
(potentially) the merger of two neutron st ars II 101 • Thus, they may trace both young and 



old stellar populations. See Refs Hi 1L I12L 1131 . 11411 for more on the development of this 
model. 

Here, we briefly review the millisecond proto-magnetar model for long-duration 
GRBs. We then discuss some of the particulars of the model, including some of its lim- 
itations and uncertainties. A recent comparison between the collapsar and millisecond- 
magnetar models can be found in Ref H- 

THE BIRTH OF MILLISECOND PROTO-MAGNETARS 

Core-collapse supernovae (SNe) leave behind hot proto-neutron stars that cool on the 
Kelvin-Helmholtz timescale (?kh ~ 10 — 100 s) by radiating their gravitational binding 
energy (~ io 53 - 5 ergs) in neutrinos |fl6l \vft\ . A fraction of these neutrinos deposit 
their energy in the tenuous and extended atmosphere of the proto-neutron starQ In the 
standard picture, net neutrino heating drives a thermal wind that emerges into the post- 
SN shock environment, blowing a wind-driven bubble into the exploding and expanding 
SN cavity [|18L ll9ll. A neutrino-driven wind is obtained in all successful models of SNe, 
regardless of how the explosion is seeded and is a natural consequence of the low- 
pressure cavity created by the explosion and the very high thermal pressure at the proto- 
neutron star surface |fl9l I20L I2H1FI For most massive stellar progenitors with extended 
hydrogen envelopes (Type-II), the cooling phase is over well before shock breakout at 
the progenitor's surface (~ 1 hour after collapse). For compact Type-Ibc SNe, the SN 
Shockwave traverses the progenitor on a timescale comparable to ?kh- 

For a neutron star with M = 1.4M , R = 10 km, rotation frequency Q, and a 
magnetar- strength magnetic field, the mass-loss rate during ?kh is given by [[22, 23, 



M(t) «4x 10- 7 4^ L5 (0exp[ft(0 2 /«c] M © s ^ (D 

where Ly.Sl.sM = L v (t) /10 51 - 5 ergs s _1 is the total neutrino luminosity, and L v « (fiv) 4 
has been assumed, where (e v ) is the average neutrino energy. The normalization of 
equation (OQ) and its L v dependence follow from the physics of the weak interaction and 
the depth of the NS gravitational potential 1I22L I23I1. The exponential factor in equation 
(OQ) accounts for the mass loss enhancement by magneto-centrifugal forces when B and 
Q are large (H c « 2000L° ° 5 8 1 5 rad s" 1 [0,H). 

One of the most important components of the proto-magnetar model is that as the neu- 
tron star cools, L v (t) and (e v ) (t) decrease on a timescale ?kh- Typically, the dependence 
L v (t) t~ l is found in cooling models of non-rotating non-magnetic proto-neutron stars 
until t ~ 30 — 40 s when L v (?) plummets as the neutron star becomes transparent tl7|] . As 
a result of equation (OQ), as L v (t) decreases, M(t) decreases. For this reason, for fixed sur- 
face magnetic field strength B, the wind becomes increasingly magnetically-dominated 



The charged-current interactions v e n -H- pe and v e p -H- ne + dominate heating and cooling. 
This statement is true even if the wind is "one-sided" in its first seconds, as in Ref HH. 



and relativistic as a function of time. Like beads on a stiff wire, the matter is accelerated 
off of the surface of the proto-neutron star by the combined action of big B and Q.. The 
degree to which the magnetic field dominates the dynamics and accelerates the flow is 
quantified by the magnetization at the light cylinder (Ri = c/Q, ~ 50Pi km): 

<r LC = B 2 /(Anpc 2 ) | LC « 15B\ 5 P^L-%{t), (2) 

where Bi$ = 5/10 15 G is the surface dipole field strength. The Lorentz factor of an 



outflow at large distances can approach — > Olc (see Ref Q24Q ) if there is efficient 
conversion of magnetic energy into kinetic energy (in which case a(r) itself will become 
small at large r). Thus, under the assumption of efficient conversion of magnetic energy 
to kinetic energy y^{t) increases dramatically as M drops on ?kh- 

Starting from a time t ~ 1 s after the collapse and successful SN explosion, there are 
several phases in the life a proto-neutron star wind. Although the specific timing of each 
phase depends on B and Q., for a millisecond magnetar they are roughly the following. 
For a second or so after explosion, as the proto-neutron star contracts to its final radius 
and spins up, the wind is pressure-dominated and driven by neutrino energy deposition. 
The asymptotic velocity of the flow is ~ 0.1c. Just a few seconds later, the wind 
rapidly becomes magneto-centrifugally dominated, but it is still non-relativistic (etc < 
1). Perhaps ~ 2 — 5 s into the cooling epoch, M(t) decreases sufficiently that Olc(?) 
becomes larger than unity, the flow is accelerated to near c by magneto-centrifugal forces 
from the surface of the neutron star out to the light cylinder. It is this late-time outflow 
— deep into the cooling epoch of the magnetar, many seconds after the explosion is 
initiated — that is the most promising for producing a GRB. Eventually, at the end of 
the cooling epoch L v (t) and M(t) decrease dramatically (see Ref flJ/ZD), o(t) increases to 
~ 10 6 and the magneto-hydrodynamical mass loss ceases, and the millisecond magnetar 
enters its "pulsar" phase in which its energy loss rate is presumed to be given by the 
force-free ("vacuum dipole") spindown formula [2511 

E FF ^B 2 R 6 a 4 /c 3 . (3) 

The wind does not emerge into vacuum. It follows the preceding, SN Shockwave, 
which propagates through the progenitor at ~ 10,000 km s . As the wind becomes 
increasingly fast a few seconds after the explosion is initiated, it shocks on the slower 
outgoing SN Shockwave. This interaction produces a wind-driven "magnetar wind neb- 
ula," and the residual azimuthal component of the magnetic field within the cavity exerts 
a "pinch" force that leads to a prominent and strong jet that pun ches through the star 
||26L I27L l28|l . This effect is equivalent to that discussed in Refs [|29L 130(1 • See also H. 



UNCERTAINTIES, LIMITATIONS, & SOME PARTICULARS 

Efficiency & Internal Shocks 

As discussed above, Olc(0 an ^ (potentially) y^t) increases monotonically as a 
function of time in the proto-magnetar model. Although there may be eruptions from 
the magnetar surface during cooling which rapidly modulate M and o, the overall trend 



from low o (< 1) to high o (> 10 5 ) on a ~ 10— 100 s timescale is unavoidable. If 
o(t) Y<=o(t), this leads to very high radiative efficiency via internal shocks (app roaching 
~ 0.3; see Ref [QlJ for discussion in the context of 060614-like GRBs l3lL l32IO. 



Jet Formation 

Perhaps counter-intuitively, the models of Refs [26|, |27|,|28] show that the otherwise 



predominantly equatorial proto-magnetar wind can be efficiently converted into a jet. 
Importantly, these works find that little of spindown (wind) power is transferred to the 
"spherical" SN component of the explosion. While the slow SN shock acts to contain 
the wind and allow the jet to develop, a variety of simulations with different parameters 
(B and Q.) show that the majority of the wind power escapes along the axis of these 
axisymmetric relativistic MHD simulations. This property of the system preserves the 
overall degree of relativity of the outflow from the light cylinder to the edge of the 
progenitor (no spindown energy is "lost" to the quasi- spherical SN). It is then difficult 
to argue that the relativistic magnetar outflow acts to super-energize the SN explosion 
for comparison with events like SN 2003dh (GRB 030329; 0311 ). which may require Ni 
yields and ejecta velocities higher than typical Ic SNe ft2$ \. 

Several theoretical questions remain: (1) What happens to the process of jet formation 
in 3D, when the magnetic and spin axes of the proto-magnetar may be misaligned? (2) 
If jet formation is indeed robust, is the alternating field ("striped" wind) configuration 
expected for an oblique rotator preserved along the pole in the jet? That is, is the jet 
stable? This issue bears critically on the efficacy of magnetic reconnection models for 
the jet acceleration and emission $2$\ . 

The Maximum Energy 

The gravitational binding energy of a neutron star is of order E gmw ~ GM 2 /R ~ 
5 x 10 53 ergs. For maximal rotation, one might guess that the total amount of rota- 
tional energy that can be stored in a proto-neutron star at birth is ~ E giav , and thus 
that the maximum possible energy for a GRB powered by the rotational energy of a 
neutron star is E rot ~ E gmv . However, at very rapid rotation rates we expect strong de- 
formations of the neutron star (as in, e.g., Ref [34]), which will emit copious grav- 
itational waves. These losses will effectively spin down the proto-magnetar. For ex- 
ample, taking the kinetic power in the outflow to be £wi n d ~ MR 2 Q. 2 /2 as the flow 
starts to become magnetically-dominated, and taking the energy loss rate in gravity 
waves to be £ G w = (32 / '5) GI 2 £ 2 £l 6 /c 5 Wl, where / = 7 45 1 45 cgs is the moment of 
inertia and e = Co.oiO.01 is the ellipticity, the critical Q. above which £gw > £wmd is 

Q ~ 500Qafi^}^i^ 1/2 %oi 2rad s_1 ' where ^-3 = M/IO- 3 M© s" 1 and R l0km = 
i?/10km. Note that the magnetic field within the neutron star may produce non-zero 
e [3^,[37]. Although sophisticated models of rapidly rotating neutron star birth are re- 



quired to address this question in detail, this estimate and others suggests that the mag- 
netar can have a maximum total rotational energy of E max ~ 5 x 10 52 ergs at the time 
of cooling. Such a calculation is particularly relevant since recent work suggests that 



some GRBs may approach this bound N38L I39TI . However, the inference of the total en 



ergy of the GRB depends on observation of the jet "break" in the afterglow lightcurve. 



Work by Ref [40] suggests earlier calculations of the total energetics via jet breaks may 
overestimate the energy in the GRB jet by a factor of ~ 3 — 4. 

A theoretical investigation of rapidly rotating 3D MHD stellar collapse is required to 
address the issue of E max for millisecond proto-magnetars in detail. Although one might 
guess that there is no similar bound for the collapsar mechanism (short of the entire 
rest mass energy of the accreted material), Ref |6(] have found that it might be difficult 
to promptly form very rapidly rotating black holes in stellar collapse because of rapid 
mass loss. Nevertheless, E max remains an important potential observational probe of the 
GRB central engine. It is particularly interesting to speculate that observations may one 
day provide evidence for a well-defined maximum energy for long-duration GRBs, thus 
providing a crucial clue to their nature. 



The Maximum Duration 

Equation © shows that Clc (and, hence, y^) is a strong function of L v . In turn, L v is 
a strong function of time in non-rotating models of proto-neutron star cooling ||l6L ll7|l . 
This time evolution implies that the wind quickly evolves from o ~ 1 to o ~ 1000 on a 
timescale of less than ~ 30 s. If the prompt emission mechanism requires that o is less 
than some upper bound (say, 1000; see 114 ill ), then the sharp cutoff in L v at the end of the 
cooling phase should quickly end the prompt emission. This implies that proto-magnetar 
powered prompt emission from GRBs should have a sharp temporal distribution in the 
rest frame (see |42|]). 

In the magnetar model, can the duration of the prompt emission be increased? There 
are potentially two ways. First, the neutrino diffusion timescale is longer for more 
massive NSs. For a 2M NS, the time at which L v suddenly decreases, increases 
somewhat, from ~ 30 s for a 1 .4 to perhaps 40 — 50 s. The second way is through rapid 
rotation. Although models of rapidly rotating neutron star cooling have not yet been 
constructed, parameterized ID models of core-collapse SNe with rotation indicate that 
the overall neutrino luminosities and average energies might be decreased by a factor of 
at most ~ 5 — 6 [14311 . Since, at zeroth order, the same amount of gravitational binding 
energy must be liberated, we expect that ?kh should increase by the same factor. Thus, 
if an upper bound on a (or y^) determines when the prompt emission shuts off, then 
the maximum duration of a long-duration GRB powered by a ms-magnetar is no longer 
than ~ 200 s. Observations of GRBs showing that the prompt phase lasts more than this 
in the rest frame would strain the magnetar model. 



Flares & Late-Time Activity 

At late times, after neutrino-driven mass loss from the proto-neutron star has subsided 
and deleptonization is complete, the energy loss rate of the magnetar should approach 
the force-free limit, given in equation ©. This implies that E Q. 4 oc t , However, 
we know from studies of pulsars in the Galaxy that the observed braking index is 
n = Q.Q./& 2 ~ 2.5 — 2.9 (see BSD, and references therein). While the cause of n ^ 3 
braking index is unknown, this finding implies that we should expect expect E oc t~ 233 
to oc r~ 2 05 in the force-free pulsar-like phase, if newly born magnetars spin down like 
their less-energetic pulsar cousins. In any case, the rotational energy reservoir should 



likely decrease monotonically, and this behavior should imprint itself on any emission 
during this phase. For this reason, the late-time flares have been analyzed by Ref 114411 and 
ll45ll . The former find that the average flare luminosity declines as a powerlaw in time 
as L f-i- 5±016 5 whereas the latter find that the peak luminosity of the flares correlates 
with the timescale of the peak asL« ? -2.7±0.5 f ormer [ s likely too shallow a time- 
dependence for the magnetar model if the flares are interpreted as powered by spindown 
||44tl . whereas the latter may be somewhat too steep. 



Recent work by Ref 114611 and 114711 explicitly connect the X-ray plateaus seen in the 
early afterglow of many Swift long-duration GRBs with magnetar spindown. If this 
correspondence is correct, the magnetar model should predict correlations between the 
character of the prompt emission and the duration and flux of the plateau. 

Also regarding very late-time central engine activity, it has been suggested that the 
discovery of a giant magnetar flare (a scaled-up version of those that occur from Galactic 
magnetars) associated with an old GRB afterglow in a nearby galaxy would provide a 
striking confirmation of the magnetar GRB model ||4~8|l . 

Time Evolution of the Magnetic Field 

A complicating and uncertain feature of the magnetar model for GRBs is the time 
evolution of the magnetic field. This could be important for both the early-time and 
late-time evolution and emission. At early times, during the neutrino cooling phase, B 
may evolve because of dynamo action in the_proto-magnetar during the explosion itself, 
or during its convective cooling phase [4^, [50]. This effect could extend or shorten the 



prompt phase if a certain range of Olc ( or Y°°) * s required for efficient emission. After 
the prompt phase, B might evolve in the magnetar via magnetic dissipation |5 hi . which 
could potentially lead to extended late-time X-ray activity. A time- variable B would also 
complicate the time dependence of spindown at all phases (eq.[3]). 

(Non-)Association with Supernovae 

Some proto-magnetars may be formed from the collapse of a Type-Ibc progenitor. In 
this scenario, one imagines that the SN mechanism in some form (e.g., the "neutrino 
mechanism" [52]) launches a SN shock with ~ 10 51 ergs that produces 56 Ni via explo- 
sive nucleosynthesis. This would explain the presence of an accompanying SN, but not 
the bright character of some GRB-SNe. If 56 Ni in excess of that produced in non-GRB 
Type-Ibc SNe is indeed required in some cases (see (530) there are two ways in which 
it might be produced by proto-magnetars: (1) some of the initial rotational energy of the 
core may be tapped rapidly via magnetic stresses, enhancing the SN shock energy as it 
is launched II43L F54TI and/or (2) as the initial slow SN Shockwave is moving outward it is 
shocked by the subsequent, highly-energetic proto-magnetar wind, again enhancing the 
shock energy [7]. Depending on the angular distribution of the wind kinetic energy, the 
latter option requires that > 10 51 ergs is extracted from the proto-magnetar on a 1 — 2 s 
timescale (see Ref 0,8). 



Proto-magnetars (and their GRBs) may also be formed by the accretion-induced 
collapse (AIC) of white dwarfs and (perhaps) by the merger of white dwarfs. In this 
scenario there is no explosive nucleosynthesis, and perhaps little 56 Ni yield JsLQli. A 



standard SN is not expected to accompany AIC due to the low total ejecta mass and 
modest quantity of radioactive ejecta j|55ll . However, Ref ll56ll show that the yield of 
56 Ni in AIC is significantly enhanced in the case of rapid rotation due to winds from 
the accretion disk that forms around the neutron star. The radioactive decay of the 
ejecta produces an optical transient that is somewhat dimmer and evolves much faster 
(on a timescale ~ 1 day) than a normal SN [[561 but which could in principle be 
detected with rapid and deep follow-up of a nearby event. Ref [58] predict that magnetar 
formation via AIC ejects up to ~ 0. 1A/ Q in highly neutron-rich material. The heavy r- 
process elements produced in such ejecta contribute a comparable amount of radioactive 
heating to 56 Ni on timescales ~ 1 day ll59ll and thus may also contribute to SN-like 
emission following the GRB. 

Importantly, if proto-magnetar-driven GRBs can arise from AIC, then long-duration 
GRBs could trace both young and old stellar populations. This is qualitatively different 
from the collapsar model. 



Does Magnetar Birth Always Produce a GRB? 



No. Estimates suggest that ~ 10% of all core-collapse SNe produce magnetars [60] 



Thus, the magnetar birth rate likely significantly exceeds the beaming-corrected GRB 



rate [610. There are at least four possibilities for relieving this potential tension: (1) many 
magnetars are born slowly rotating and thus do not produce GRBs, (2) many millisecond 
magnetars are born in Type-II progenitors with extended hydrogen envelopes and their 
jets are quenched ||62|l . (3) many millisecond magnetars simply do not form jets able 
to escape the star (because of, e.g., instabilities), and finally (4) (a variant on the first), 
because E (e.g., eq. [3]) is a very strong function of both B and Q,, only the magnetars 
born at the extremes of these distributions (e.g., B ^ 10 15,5 G, P ^ 2 ms) produce bright 
cosmological GRBs. The rest give rise to less energetic events like XRFs or normal SNe. 
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